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For many practical spacecraft applications, algorithms for determining spacecraft attitude must combine 
inputs from diverse sensors and provide redundancy in the event of sensor failure. A Kalman filter is 
suitable for this task, however, it may impose a computational burden which may be avoided by sub 
optimal methods. A sub optimal estimater is presented which was implemented successfully on the Delta 
Star spacecraft which performed a 9 month SDI flight experiment in 1989. This design sought to minimize 
algorithm complexity to accommodate the limitations of an 8K guidance computer. The algorithm used is 
interpreted in the framework of Kalman filtering and a derivation is given for the computation. 


INTRODUCTION 

Historically, satellite attitude determination has 
relied on simple deterministic calculations for 
batch processing of telemetry data because real- 
time recursive algorithms such as Kalman filters 
imposed an impractical computational burden. 
This burden has become less daunting with 
advances in flight-qualified microprocessors, 
however, simple algorithms remain important for 
maintaining the reliability and controlling the 
development cost of real-time software. 

This paper examines the algorithm used to 
estimate attitude for Delta Star. This algorithm 
applies deterministic gains to measurement 
data. Nonetheless, it is desirable to perform an 
statistical error analysis. The attitude estimation 
problem is cast as a Kalman filtering problem 
such that the performance of the sub optimal 


deterministic gains can be quantified. As a 
convenient byproduct, the Kalman gains implicit 
in this setup provide an alternative estimation 
procedure with only a modest increase in 
computations. 

DELTA STAR BACKGROUND 

The SDIO sponsored Delta Star spacecraft 
operated on-orbit for nine months during 1989. 
Its objectives included multi-spectral observation 
of low earth orbital phenomena against various 
earth and space backgrounds. Numerous 
pointing and tracking guidance modes required 
modest but reliable, knowledge of spacecraft 
attitude. 

The spacecraft consisted of two sections: a 
guidance and propulsion section and a sensor 
section. Each was controlled by separate 
processors designated guidance computer (GC) 
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and flight processor (FP) respectively. All 
primary GN&C functions resided in the GC. 
Because the guidance and propulsion section 
was based on the Delta launch vehicle second 
stage, the GC was a Delco Magic 352 guidance 
computer, featuring 8K of random access 
memory to accommodate data, program 
instructions and a resident operating system. In 
addition to GN&C functions, the GC flight 
program sequenced discretes to control avionics 
subsystems, processed telemetry and uplinks 
from the ground, and provided a protocol for 
communication with the FP. The limited memory 
budgeted for attitude determination made a 
simple design imperative. 





Figure 1. Delta Star Spacecraft 

QUATERNION CONVENTIONS 

In this section notation and conventions are 
developed for the quaternion q. The two 

primary coordinate frames of interest in this 
report are an inertial reference frame I, and a 
spacecraft body-fixed frame B. A coordinate 
frame is given by a triad of orthonormal basis 


vectors which obey the right-hand rule. A 
change of basis is specified by a rotation or 
direction cosine matrix Tf defined by 

y.=Tfy, (1.) 


where 


y. e * 3 
y, e * 3 


are the same vector expressed in I and B 
coordinates. 


The rotation matrix Tf can be represented by a 
quaternion q. the quaternion is a globally 
nonsingular mapping of the rotation matrix. The 
set of attitude quaternions is defined as 

Q = [(q. ,q. ) e % x <R 3 :<7 3 + |g r f = X q. > o} 

) 

where the first condition is the unit quaternion 
normality constraint and the second is a 
convention to eliminate the ambiguity of sign 
which arises because i,q»,qv) and (- q»,-q *) 

represent the same attitude. With these 
conventions, Tf can be computed from the 
quaternion q by the formula 

Tf = l + 2al~2q,Q q . (3.) 


where for a e 9? 3 , is defined as 
Q,:0t 3 x9t 

C = n.b (4.) 

c = axb 

for a.toeSR 3 . Quaternion multiplication is 
defined as follows: 
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• FIGURE 2. Delta Star Attitude Filter 


c = ab, a,beQ 

Ct = 3ab a—[av ,bv] /g i 

Cv = 3s bv + baBv + By xbv 

C = (C,,Cv)eQ 

Quaternion multiplication is important because it 
corresponds to compositions of quaternion 
rotations. That is, for coordinate frames A,B,C, if 
qi rotates A into B and <72 rotates B into C, then 
<73 = q-[q 2 rotates A into C. 

The quaternion q. representing T* , evolves in 
time according to the equation 

q = ^q(0, co) (6.) 

where to is the instantaneous angular velocity of 
the B-frame with respect to the l-frame specified 
in B coordinates. 


The inverse of q is denoted q * which is also 
called the conjugate of q [3] and is defined by 

q* = (q,,-q) , and qq* = q'q = (1,0) 

) 

Finally, a vector vie9t 3 in I coordinates is 
transformed into v B e 9t 3 in B coordinates by 

(0,v B ) = q'(0,vi)q = (0,T?v,) (8.) 

DELTA STAR ATTITUDE DETERMINATION 
The attitude sensors on Delta Star spacecraft 
consisted of five sun sensors and a dual conical 
scan horizon sensor. The five sun sensors were 
configured to provide omni-directional sun 
coverage. The horizon sensor had a 26°x26° 
field of view. These sensors provided attitude 
measurements for comparison against on-line 


219 











ephemeris. The FP edited sensor data for wild 
points and compensated horizon sensor data for 
earth oblateness before passing unit vectors for 
sun and nadir across a communication interface. 
Modern earth sensors are equipped to provide 
such compensation using embedded 
microprocessors. The Delta Star Attitude Filter 
(DSAF)design is shown in Figure 2. 

Traditional deterministic methods which compute 
a direction cosine matrix from a pair of 
independent measured vectors, such as the 
TRIAD algorithm [2], offer extreme simplicity but 
suffer from several deficiencies, 

1 . Only the current vector pair 
factor into the attitude estimate (i.e. 
noisy measurements are not averaged). 

2. Nearly collinear measured 
vectors produce dubious solutions 

3. The two vectors of a pair must 
be synchronous for a solution and cause 
complications if they arrive 
asynchronously. 

4. Measurements from different 
sensors cannot be weighted to reflect 
relative noise levels. 

A Kalman filter will eliminate these deficiencies. 
However, the computations required by such a 
filter were considered prohibitive for the Delta 
Star application. The design shown in Figure 2. 
also eliminates these deficiencies, but without 
the matrix computations required by the Kalman 
filter to propagate a covariance matrix and 
compute a gain as a function of the covariance. 


The constants a and P are design parameters 
used to control noise rejection and to weight 
measurements from sun and horizon sensors 
with respect to each other.A method for 
preforming a statistical error analysis of this 
design is presented below. Sub optimal gains 
are derived in terms of a and p. A statistical 
interpretation of these gains is given which 
provides considerations for selecting a and p. 

In Figure 2., a running estimate of attitude is 
maintained by integrating angular rates from 
gyros according to (6.). This running estimate 
denoted by q differs from q as a result of gyro 

drift and initial condition errors. 

A sun sensor produces two measurements from 
which the sun vector in B-coordinates can be 
derived. An earth nadir vector is similarly 
derived from the outputs of an horizon sensor. 
Specification of this processing will not be given 
here. These computations were performed by 
the FP for Delta Star are not formally considered 
a part of the DSAF design. 

We will distinguish between observation vectors 

and measurements of these vectors. An 
observation vector will be denoted by e * 6 ^ 

for time tk where A is a tag denoting the type of 
observation (S:sun,N:nadir}.and C denotes the 
coordinate frame in which the vector is 
expressed. 

A measurement z * k e gj3 of e A js derived 
from sun and horizon sensors. For our 
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purposes, the former is a "noisy" version of the 
latter. 

The vectors and * s are available from an 

e l k e i k 

on board ephemeris calculator. From these 
reference vectors estimates of e A denoted by 

v Bk 

are computed by 

(0,e* k ) = ?(0 ,ef k )q < 9 > 

The discrepancy between q and q is then 
estimated from the discrepancy between $a 

° k 

and z *. In Figure 2., we note that the 
discrepancy between and z *. is captured 

in the form of the cross product of these two 
vectors which is used to compute a corrective 
rate by which to improve the estimate q . 

The discrepancy between q and q will be 
defined by 


v£=Htdq; k +r\t;Cov(r\t) = Rf ( 1r > 
To begin, v * is defined by 


where is considered gaussian white noise 
with covariance matrix r* Given q, we can 
define /j£;9l 3 -> 9l 3 from 0 1 •) by 

v* = h*(8q V k)+T\k (13) 

The Jacobian matrix of h* g$ 3 _> g{ 3 is given by 

H k =Dh k (bq v k )- -2T?nl ( 1 4 ) 

where n is the derivative operator, f is 
computed from (2.) and e = ei 

B 

Then 


q = 5 qq ( 10 ) 

Because can be computed from g q v using 
the normality constraint in (2.), 5 q v will be used 
to define the attitude error. In the statistical error 
analysis below, we investigate the behavior of 
P = Cov(Sq„) for the DSAF 9 iven specified 
statistical assumptions. 

In the following, we derive a measurement 
sequence and matrix sequences 

{4>* ,Hk ,K*}*_ 1 such that under specified 
conditions the following error propagation and 
update equations apply for the DSAF: 


v* - h k ( 0 ) = H* 8 q* k +i\f+ oflSg, . \) 

) 


so that to first order 


v*=Hfa„ t +rrf (16) 

This is the linearized observation equation. 

By (5.) and (9.) if e = o over the interval 

[ tk.tu +i; 

8*7* ,♦ , =8q^, (17.) 

= , ;<bk = / 
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This is the linearized state provides and 
gives the linearized state equation Finally, 

q = lq(0,a + e); q(t k ) = qk ; q k+ 1 = q(t k ) 

(18.) 

To first order this is equivalent to the Extended 
Kalman Filter (EKF) update scheme defined by 
the procedure 


q k = q k 

r k =(*'lT' B z)q- k 

3 = lq(0,(o + e); q(t k ) = q{ ; q k+ , = q(t k ) 


where 


) 


<7* +1 = q* + i +o(x); x = f*+i-f* 

( 20 .) 

Then to summarize , (16.), (17.) and (19.) provide 
an a set of propagation and update equations in 
the same form as (1 1 .) with 

K k =cn Tg k (22 ) 

KZ=frT' Bk SS T ;S = e s Bk 

This signal generation model enables us to 
analyze the behavior of the covariance 

P k =Cov(5q vk ) ■ (23) 

The covariance for this setup propagates 
according to 

P k ^=F*i P k (F*) T + G£ 

) 

where 

F{ -{l-KtHt);Gt=KtRt(Kf) T 
By (17.) we need not distinguish between pre- 
update and post-update covarances (i.e.. 


p- - p+). If pA is constant, which we shall 
assume, then it is easy to see thang>> is also 
constant. 



HS Processing 


HS Field of 
View 


FIGURE 3. Solar Inertial Geometry 


To understand the significance of a, and P in 
(22.), consider the simple geometry shown in 
Figure 3.. The sun vector lies in the orbit plane 
and intervals of sun sensor and horizon sensor 
usage are as shown. Define the set of basis 
vectors s-BvB 2 wh e r e s is the sun vector and 
BvB 2 are chosen to form a right handed 
orthonormal triad or coordinate frame. We will 
call this coordinate from the I'-frame. 


The vectors s-BvB 2 are al1 eigenvectors of 
F f . with eigenvalues x s ,\ Bii \ B2 such that 

Xs- 1 (25.) 

Xb, = Xb 2 = 1 - ax 

The matrix pS modifies the covariance 

f k 

according to (24.) when sun sensor data is 
processed. If 5 q v is expressed in I' 

coordinates, then by (24.) 


a i* + i = 4<4 * + i + Ys 

(j2 B^ k ^=^ 2 B^B, k+ : + yB, 
CT fi2* + 1 = ^S 2 a B2* +1 +702 


(26.) 
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where 


a| k = P1.1* ; ys = Gi,i* 

a| )fc = P2.2* ; YBt = G2.2* 

a| Jk = Pa,3* ; Yfls = 63,3* 

and it is assumed that q * is constant. The 
behavior of the variances defined in (26.) is 
simple to understand in terms of the difference 
equations. The error around the sun vector 
increases at a constant rate (in terms of 
variance) at a rate determined by the sensor 
noise and the gain a. The orthogonal 
components decay to a steady-state value as 
the corresponding eigenvalues are less than 
unity. The steady-state residual can be 
computed using the Final Value Theorem for Z- 
transforms. The decay rate is exponential and 
easily determined from ,Xs 2 • A d esi 9 n value 
for a is achieved by establishing acceptable 
values for error growth around the sun vector, 
and steady-state residual and decay rate for 
error about the orthogonal vectors, and trading 
off one for the other for an “optimal" compromise 
The horizon sensor gain P, can be selected 
similarly. 


ERROR ANALYSIS 


The covariance propagation above is limited 
because only the effects of sensor noise are 
considered. To investigate the effects of other 
errors such as a constant gyro drift, the method 
described in [1] is used. The basic idea is shown 
in Figure 4. 



FIGURE 4. Error Analysis Method 


For the geometry shown in Figure 3 this method 
was used to generate a covariance history 
assuming that, 

/?* =D/afltf3.05x10 -6 ,3.05x10“ 6 ,3.05x1(T 6 ; 
f?" = Diag[8. 46 x 1 (T* , 8. 46 x 1 (T 8 , 8. 46 x 1 (T* ] 
pf = Dlag[7. 61 X 1 0 -7 , 7. 61 x 1 0" 7 , 7. 61 x 1 0~ 7 ) 

Qk = Dlag[6. 53 x 1 0" 13 , 6. 53 x 1 0 -13 , 6. 53 x 1 0" 1 3 7 
a = 0.01 
P = 0.06 

where Q k is the covariance of the constant gyro 
bias error in radians per second, quaternion 
error is dimensionless but, approximately half of 
angular error in radians and sensor error is 
similarly approximately half of the angle error 
produced by sensor noise given in radians The 
factor of two comes from the definition of 
quaternion in terms of rotation angle and rotation 
vector{3]. The result is shown in Figure 5. 
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For comparison, the gains in (22.) are replaced 
by Kalman gains computed using 

Kk = PxHf(H£P k H* T + R*J' 

(27.) 

The result is shown in Figure 6 

In each of these cases, the spacecraft begins by 
processing sun sensor data. The error around 
the sun line slowly increases, and the orthogonal 
components are reduced. At approximately 
1100 seconds, an horizon sensor update occurs. 
On the time interval (1200-4000), no sensor data 
is processed and pure gyro drift is observed. At 
4000 seconds sun sensor data is processed 
again and a new cycle begins. We observe that 
the convergence rates are faster for the Kalman 
filter, that the Kalman filter variances converge 
to smaller values and that 
the 
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FIGURE 6. Kalman Filter Covariance Analysis 

Kalman filter performs better with respect to 
orthogonal components of error during horizon 
sensor updating. Note that no attempt is made 
to estimate gyro drift from the sensor data. 



FIGURE 7. Nadir Pointing Geometry 


We see that for the scenario described by Figure 
4. the DSAF compares favorably with the 
Kalman filter without having to propagate a 
covariance or compute a Kalman gain. The 
Kalman filter does, however, afford an 
advantage which is not evident in the above 
analysis. The DSAF will not worfc if only horizon 
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sensor data is available. Clearly, such a 
capability is desirable for attitude determination 
reliability in the event of a sensor failure. The 
Kalman filter does have this capability. For the 
same statistical assumptions as above, but 
using only horizon sensor data for the geometry 
shown in Figure 7. we obtain the covariance 

history shown in Figure fl. 
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presented, and performance is compared 
against a Kalman filter. 
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FIGURE 8. DSAF Covariance Analysis 
Horizon Sensor Only 


CONCLUSION 


This paper has presented a simple filtering 
algorithm DSAF for determining spacecraft 
attitude from vector observations. This algorithm 
was used successfully on-orbit for the Delta Star 
SDIO flight experiment in 1989. It offers several 
advantages over simple deterministic methods 
such as TRIAD, but does not require as much 
computation as a Kalman filter mechanization. If 
a Kalman filter is required or desired for an 
application, the DSAF is easily extendible to a 
Kalman filter by means of a more elaborate gain 
computation. The design parameters of the 
DSAF are motivated, an error analysis is 
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